Temperature Chaos and Bond Chaos in the Edwards- Anderson Ising Spin Glasses 

Domain- Wall Free-Energy Measurements 
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Domain-wall free-energy 8F, entropy SS, and the correlation function, Ctemp, of 5F are measured 
independently in the four-dimensional ±J Edwards-Anderson (EA) Ising spin glass. The stiffness 
exponent 0, the fractal dimension of domain walls d B and the chaos exponent £ are extracted from 
the finite-size scaling analysis of SF, SS and Ctemp respectively well inside the spin-glass phase. The 
three exponents are confirmed to satisfy the scaling relation £ = d s /2 — 9 derived by the droplet 
theory within our numerical accuracy. We also study bond chaos induced by random variation 
of bonds, and find that the bond and temperature perturbations yield the universal chaos effects 
described by a common scaling function and the chaos exponent. These results strongly support the 
appropriateness of the droplet theory for the description of chaos effect in the EA Ising spin glasses. 
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PACS numbers: 75.10.Nr, 75.40.Mg, 05.10.Ln 

In randomly frustrated systems such as spin glasses, 
directed polymer in random media (DPRM) and vortex 
glasses, the equilibrium ordered state could be completely 
reorganized by an infinitesimally small change in envi- 
ronment 0, 0, 0, IS • This curious property called 
chaos effect has attracted much attention since it was 
found in 1980s Especially, chaos induced by tem- 
perature variation (temperature chaos) is now of great 
interest because of its potential relevance for rejuvena- 
tion caused by temperature variation 0|. However, the 
issue of temperature chaos still remains far from be- 
ing resolved. In particular, concerning low-dimensional 
Edwards- Anderson (EA) Ising spin glass models, the sit- 
uation is very controversial because numerical studies so 
far done provide the evidence both for and against tem- 
perature chaos 0,0, HE Ell- 
in the present work, we examine temperature chaos by 
numerical measurements of the domain wall free-energy 
SF, the difference in the free-energy between the system 
with the periodic boundary condition (BC) and that with 
the anti-periodic BC. This SF relates to the effective cou- 
pling Jeff between the two boundary spins Sl and Sr (see 
Fig.[TJ as J c ff = -5F/2 [HO]!. We find 6F of each sam- 
ple exhibits oscillations along the temperature axis pro- 
viding direct evidence of the temperature chaos. Further- 
more, we find from simultaneous observations of the do- 
main wall energy SE and so entropy SS that SE and TSS 
are large but they cancel with each other in the leading 
order to yield significantly small SF = SE — TSS. Such 
intriguing behavior is indeed predicted by the droplet 
theory [14|. For a quantitative check of the droplet the- 
ory we focus on the anticipated scaling relation 



derived from it, where the stiffness exponent 9 is ex- 
tracted from of ~ L e , the fractal dimension of domain 
walls d s from as ~ L da / 2 , and the so called chaos expo- 
nent C from Ct cm p (L, T,T + A) ~ /(L(AT)- 1 /?). Here 
uf and as are the standard deviations of SF and SS, re- 
spectively, and Ctemp is the correlation function of SF's 
defined by Eq.(@J below and f(x) is a certain scaling 
function. We find the three fundamental exponents thus 
extracted indeed satisfy Eq. 0) well inside the spin-glass 
phase whose thermodynamic properties are dominantly 
governed by the T = fixed point. 

We also study bond chaos by measuring how SF varies 
with changes in couplings. The result evidently shows 
the existence of bond chaos. Moreover, the scaling anal- 
ysis of two correlation functions associated with temper- 
ature and bond perturbations reveals quantitatively that 
not only the chaos exponent but also the scaling func- 
tion are common to both the perturbations. This uni- 
versal aspect of chaos effect anticipated from the droplet 
theory Mjl is also observed in the Migdal-Kadanoff spin 
glasses |l5l and the DPRM |(| . All of our numerical 
results, particularly the quantitative check of Eq. QJ, arc 
strong evidence not only for the existence of chaos in the 
EA Ising spin glasses but also for the appropriateness of 
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C = d s /2 - 9 
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FIG. 1: Model for the boundary flip MC. 
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the droplet theory for its description. 

The boundary flip MC method — Let us first describe 
the boundary flip MC method 0, 0] which enables us 
to measure the domain-wall free-energy. We consider a 
model which consists of Ising spins on a ci-dimcnsional 
hyper-cubic lattice of L d and two boundary Ising spins 
5l and Sr. (see Fig.^l. The usual periodic BC is applied 
for the directions along which the two boundary spins 
do not lie. The Hamiltonian is H = — ^luj) JijSiSj, 
where the sum is over all the nearest neighboring pairs 
including those consisting of one of the two boundary 
spins and a spin on the surfaces of the lattice. In our 
boundary flip MC simulation, the two boundary spins 
are also updated according to a standard MC procedure. 
For each spin configuration simulated, we regard the BC 
as periodic (anti-periodic) when 5l and Sr are in par- 
allel (anti-parallel). Since the probability P P (ap)(T) for 
finding the periodic (anti-periodic) BC is proportional to 
cxp[— F P (ap)(T)/T], where Fp(ap)(P) is the free-energy 
with the periodic (anti-periodic) BC, we obtain, with 
Pap = 1--Pp, 

SF(T) = F P (T)~F AP (T) 

= -fc B T{log[Pp(T)]-log[l-P P (T)]}. (2) 

We also measure the thermally averaged energy 
£p(AP) (T) when the two boundary spins are in parallel 
(anti-parallel). It enables us to estimate the domain- wall 
energy SE(T) = E P (T)-E AP (T). Then, the domain- wall 
entropy SS is evaluated either from SS = (5E — 8F)/T or 
SS = — d^- 1 . We have checked that both the estimations 
yield identical results within our numerical accuracy. 

We study the four-dimensional ±J Ising spin glasses 
in the present work. In four dimensions the value of 
the stiffness exponent 6 is significantly large 0, Il9j . 
which enables us to make scaling analyses rather eas- 
ily as compared in three dimensions. The values of {Jij} 
are taken from a bimodal distribution with equal weig hts 
at = ±J. We use the exchange MC method |2fJ to 
accelerate the equilibration. The temperature range we 
investigate is between 0.6 J and 4.5 J, whereas the critical 
temperature of the model is around 2.0 J |21|. The sizes 
we study are L = 4, 5, 6, 7, 8 and 10. The number 
of samples is 824 for L = 10 and 1500 for the others. 
The period for thermalization and that for measurement 
are set sufficiently (at least 5 times) larger than the er- 
godic time, which is defined by the average MC step for 
a specific replica to move from the lowest to the highest 
temperature and return to the lowest one. 

Temperature chaos — In Fig. [21 we show temperature 
dependence of 5F, SE and TSS for 5 samples. Oscil- 
lations of the three observables become stronger with 
increasing L. We in fact see that SF of some samples 
changes its sign, meaning that the favorable BC with the 
lower free-energy changes with temperature. We also see 
that, as predicted by the droplet theory 0], SE(T) and 



T5S(T) exhibit very similar temperature dependence and 
cancel with each other in the leading order to yield rela- 
tively small SF. 

In Fig. the standard deviations, op, 0e and its, at 
T = 0.6J are plotted as a function of L. Interestingly, 
<7g, which gives the amplitude of | | , increases more 
rapidly than of, i.e., the amplitude of SF. See 0] for a 
similar observation in the three-dimensional EA model. 
As argued by Banavar and Bray |l5l | , this result naturally 
leads us to the conclusion that SF in the limit L — > oo is 
totally temperature chaotic. 

The inset of Fig. shows 0{T) and d s (T)/2 estimated 
by linear least-square fits of ln(aF) and ln(as) against 
ln(L) at each temperature. As expected from the droplet 
theory which is constructed around the T = fixed point, 
the two exponents converge to a certain value at low tem- 
peratures. By averaging over the lowest five tempera- 
tures, we obtain 

6 = 0.69 ±0.03, d s = 3.42 ±0.06. (3) 

Our 9 is compatible with other estimations while 
our d s is somewhat smaller than other ones |22j. The 
apparent temperature dependence of 0(T) and d s (T) at 
higher temperatures is considered to be due to the crit- 
ical fluctuation associated with the unstable fixed point 
at T c , combined with the finite-size effect. Its detailed 
quantitative analysis is, however, beyond the scope of 
the present work. 
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FIG. 2: SF (left), SE (middle) and TSS (right) vs. tempera- 
ture for 5 samples. L = 4, 6, 8 and 10 from top to bottom. 
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We next examine the correlation function defined by 



temp 



(L, T,T + AT) 



5F{L,T)5F{L,T+AT) 
a F (L,T)a F (L,T + AT) ' 



(4) 



where 7TT is the sample average. A similar correlation 
function was first introduced by Bray and Moore to study 
bond chaos Q . The inset of Fig. 0] shows the raw data 



of 1 — Ctemp at T = 0.6 J. Ctcmp approaches zero rapidly 
with increasing L. From the prediction of the overlap 
length by the droplet theory, over which the configura- 
tions at the two temperatures are unrelated, we expect 
one parameter scaling of Ct om p = /(LAT 1 ^) whose test 
is shown in the main frame of Fig. We see that the 
scaling works nicely. The value of £ is evaluated to be 
1.12 ± 0.05 by the fitting. Quite interestingly, this value 
is consistent with the value £ = 1.02 ± 0.06 obtained by 
substituting eq. into Eq. predicted by the droplet 
theory. This is one of the main results of the present 
work. We also see that the data are consistent with the 
expected asymptotic behavior in the limit L^AT — > 0, 



1-Q 



temp 



(L^AT) 2 jj], as depicted by the line. 



Bond chaos and universality — We also study bond 
chaos by comparing two systems with correlated coupling 
sets. The perturbed couplings {J'ij} are obtained from 
the unperturbed ones {Jij} by changing the sign of J„ 
with probability p. Since simulation for bond chaos costs 
much more time than that for temperature chaos, we only 
examined L = 4, 6, 8 for bond chaos. 

Now let us consider an observable 5S hond = — SF ^f F , 
where AJ = ^Jp and 8F (SF') is the domain-wall free- 
energy of the unperturbed (perturbed) system. SS hond 
here and SS discussed above are similar in a sense that 
the both are the increment ratios of SF against the per- 
turbations. The ratio against AJ, not p itself, is con- 



sidered to compare temperature perturbation and bond 
perturbation properly [5j. In Fig. [21 the standard devia- 
tion of 5S hond , denoted as cr£ ond , is also shown. SS hond is 
estimated with AJ ps 0.03, which corresponds to a small 
value of p ps 0.0009. The line for a bond and that for a s 
have the same slope, which suggests that temperature 
and bond perturbations belong to the same universality 
class. The coefficient of crg ond is, however, about 16.4 
times as large as that of as- In the inset of Fig. 0we 
show the raw data of the correlation function for bond 
perturbation defined by 



Cbond{L,T,p) 



_ 5F(L,T)5F'{L,T) 
= a F (L,T)a' F (L,T) ' 



(5) 



Again, the correlation decays faster with increasing L. 
In the main frame of Fig. [SJ we test a similar scaling 
to that in Fig. 21 by assuming that the overlap length of 
the bond perturbation scales as AJ -1 /^. All the data 
again collapse into a single curve. The chaos exponent £ 
is evaluated to be 1.10 ± 0.10 by the fitting. 

To compare the two scaling functions for temperature 
and bond perturbations, we plot in Fig. [5] all the data of 
both Ctemp and Cbond by using the same chaos exponent. 
Here we use the value £ = 1.12 in Fig. 21 and multiply the 
scaling variable L(AJ) 1 ^ by factor 17.5. All the data 
roughly merge into a single curve, indicating that the 
chaos exponent and the scaling function for temperature 
chaos are the same as those for bond chaos. Lastly, by 
estimating the overlap length £ as the value of L for which 
C = 0.5, we obtain 



t-temp 



11.5AT 



-i/C 



0.657AJ 



-i/C 



(6) 



where £ ps 1.12. We see that the overlap length of 
the bond perturbation is much shorter than that of the 




FIG. 3: Size dependences of <tf, <te, cs and o"s° n at T = 0.6J. 
See text for their definitions. The two straight lines for of 
and as are obtained by linear least-square fits of ln(<TF)/ln(crs) 
against ln(L). The line for o"g ond has the same slope as that 
for as- The inset shows the data for f9' cff ) and d!f S ^ /2. 
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FIG. 4: A scaling plot of 1- Ctemp (i, T, T+ AT) at T = 0.6 J 
against L(AT) 1 ^ with £ = 1.12. The line is proportional to 
x 2(t . In the inset, 1 — Ctcmp for L — 4, 5, 6, 7, 8 and 10 (from 
bottom to top) are plotted as a function of AT. 
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FIG. 5: A scaling plot of 1 - C bond (L,T,p) at T = 0.6J 
against L(AJ) 1/C with ( = 1.10, where AJ = ^/p. The line 
is proportional to x 2< * . In the inset, the raw data for L = 
4, 6 and 8 (from bottom to top) are shown as a function of p. 



thermodynamic properties are dominantly governed by 
the T = fixed point. These results are certainly strong 
evidences for the appropriateness of the droplet theory 
for the description of chaos effect in the EA Ising spin 
glasses. On the other hand, recent work by Rizzo and 
Crisanti indicates the existence of similar chaos ef- 
fects in the Sherrington-Kirkpatrick model. Whether our 
results are consistent with the mean field view point or 
not is an interesting open problem. 
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FIG. 6: A scaling plot of 1 — C by using all the data in 
Figs.|l|and|K| The scaling variable for the data of temperature 
perturbation is L(AT) 1 ^ , while that for the data of bond 
perturbation is 17.5 x L(AJ) 1 ^, where £ = 1.12. The line is 
proportional to x 2C ' . The inset shows the same plot for C. 

temperature one. This result as well as factor 16.4 
in the 0\g ond scaling mentioned above are the quantita- 
tive description of the well-known fact that temperature 
chaos is much more difficult to be observed than bond 
chaos HEIl. 

Conclusion — In the present work, we have studied 
the four-dimensional EA Ising spin glass with focus on 
the chaos effect. As a consequence, many non-trivial 
predictions by the droplet theory, such as the SF oscil- 
lation along the temperature axis and the cancellation 
of 5E and TSS are found in this model. Most impor- 
tantly, the scaling relation of Eq. Q and the universal 
aspect of temperature and bond chaos effects are quanti- 
tatively confirmed well inside its spin-glass phase whose 
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